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ABSTRACT 

We report the discovery of PSR J1909— 3744, a 2.95 millisecond pulsar in 
a nearly circular 1.53 day orbit. Its narrow pulse width of 43 /is allows pulse 
arrival times to be determined with great accuracy. We have spectroscopically 
identified the companion as a moderately hot (T ^ 8500 K) white dwarf with 
strong absorption lines. Radial velocity measurements of the companion will yield 
the mass ratio of the system. Our timing data suggest the presence of Shapiro 
delay; we expect that further timing observations, combined with the mass ratio, 
will allow the first accurate determination of a millisecond pulsar mass. We have 
measured the timing parallax and proper motion for this pulsar which indicate a 
transverse velocity of 140j]4o km s~^. This pulsar's stunningly narrow pulse profile 
makes it an excellent candidate for precision timing experiments that attempt 
to detect low frequency gravitational waves from coalescing supermassive black 
hole binaries. 

Subject headings: binariesxlose — pulsars: individual(PSR J1909— 3744) — rel- 
ativity — stars: distances — stars: neutron — white dwarfs 



1. INTRODUCTION 

Binary radio pulsars provide a rich laboratory for a wide range of physical inquiry, 
including neutron star masses and the evolution of binary systems. The pulsars with the best- 
determined masses all have eccentric relativistic orbits, spin periods of tens of milliseconds, 
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and masses clustered very tightly around 1.35 Mq (Thorsett & Chakrabarty 1999; Bailes et al. 
2003). It is generally thought that millisecond pulsars (P < 10 ms) should have greater mass 
than these longer-period pulsars due to a larger accreted mass. There is indeed a statistical 
suggestion that this is the case, but to date no precise millisecond pulsar mass measurements 
exist. Currently, the best-determined millisecond pulsar masses are nip = l.ST^gJ^M© for 
PSR B1855+09 (Nice et al. 2003) and rup = 1.58 ± 0.18 for PSR J0437-4715 (van 
Straten et al. 2001). 

PSR J1909— 3744 was discovered during the Swinburne High Latitude Pulsar Survey, a 
recently-completed pulsar search using the 13-bcam multibcam receiver on the 64-m Parkcs 
radio telescope (Jacoby 2003; Jacoby et al. in preparation). This survey is similar to, and an 
extension of, the highly-successful Swinburne Intermediate Latitude Pulsar Survey (Edwards 
et al. 2001; Edwards & Bailes 2001a,b). The pulsar was discovered in three adjacent survey 
beams observed on 2001 January 25 and confirmed on 2001 May 26. Due to the pulsar's 
strong broad-band scintillation, the initial position inferred from the discovery data was in 
error by nearly a full beamwidth; it was not until we obtained a phase-connected timing 
solution that we realized the correct position. 

Among the properties that can make a pulsar useful for precision timing experiments 
are a bright, narrow pulse profile, extremely stable rotation, proximity to the earth, and a 
companion star which can be detected and studied. PSR J1909— 3744 appears to possess all 
of these desirable traits. Nearly all of the pulsar's fiux comes in a single narrow, sharp peak 
with full-width half-maximum (FWHM) of only 43 /xs, less than 1.5% of the 2.95 ms pulse 
period (Fig. 1). The brightness of the pulsar varies by at least a factor of 30 in the 21-cm 
band. The mean flux density of our observations is 3mJy, but this is biased by the fact 
that we tend to observe longer when the pulsar is bright. Extremely high precision timing 
observations are possible during episodes of favorable scintillation. 

2. PULSE TIMING 

Following the confirmation of PSR J1909— 3744, we began timing observations with 
the 2 X 512 x 0.5 MHz Parkes fllterbank centered on 1390 MHz. As a result of our incor- 
rect discovery position, few observations taken before February 2002 yielded data useful for 
high-precision timing. Integration times ranged from three to 100 minutes, with the longer 
observations broken into 10-minute sub-integrations. 

We followed standard pulsar timing procedures: folded pulse proflles from individual 
observations were cross-correlated with a high signal-to-noise template profile to determine 
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an average pulse time of arrival (TOA), corrected to UTC(NIST). The standard pulsar timing 

package tempo^, along with the Jet Propulsion Laboratory's DE405 ephemeris, was used for 
all timing analysis. All 1.4 GHz TOAs with uncertainty less than 1 //s were included. TOA 
uncertainties were multiplied by 1.25 to achieve reduced — 1- Several observations were 
also obtained in the 50-cm band with the 2 x 256 x 0.125 MHz Parkes filterbank centered 
on 660 MHz. As these observations did not yield TOA uncertainties less than 1 /is they were 
used only to determine the dispersion measure (DM), which was then held fixed. 

Because of the system's low eccentricity (e), we used the ELLl binary model which 
replaces the familiar longitude of pcriastron (u), time of pcriastron (Tq), and e with the time 
of ascending node (Tgsc) and the Laplace-Lagrange parameters esina; and ecosa; (Lange 
et al. 2001). We give the results of our timing analysis in Table 1 and show residuals relative 
to this best-fit model in Figure 2. The weighted RMS residual is only 330 ns, among the 
best obtained for any pulsar over ~2 year timespans. 

We have begun regular observations of PSR J1909— 3744 with the Caltech-Parkes- 
Swinburne Recorder II (CPSR2; Bailes 2003). CPSR2 is a baseband recorder and real-time 
coherent dedispersion system with two dual-polarization bands, each 64 MHz wide. The 
CPSR2 pulse profile shown in Figure 1 is smeared by only 2//s, compared with the 82 /xs 
time resolution of the 512-channcl filterbank. Initial results from CPSR2 data show great 
promise for high precision timing, and will be reported in a future paper. 

2.1. Shapiro Delay 

Though our current data set suffers from poor orbital phase coverage (specifically a lack 
of high-quality observations near inferior conjunction), we find strong evidence of Shapiro 
delay. If we remove Shapiro delay (i.e. companion mass, rric, and sine of orbital inclination, 
sini) from the timing model and allow TEMPO to optimize the remaining parameters, 
increases from 98.8 to 128.1 (for 102 and 104 degrees of freedom respectively), indicating a 
highly significant detection. The pcak-to-pcak residual is essentially the same in either case 
(3.4 /is) as much of the Shapiro delay signal is absorbed into an erroneously large Roemer 
delay; however, the weighted RMS residual increases to 370 ns when the model excludes 
rric and sini. Holding the astrometric, orbital, and spin parameters fixed at the values in 
Table 1 but neglecting Shapiro delay increases the peak-to-peak residual to 8.1 //s, with 
pulses observed nearest inferior conjunction arriving later than predicted by the model. 



^http: / /pulsar. princeton.edu/tempo 
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Given its large uncertainty, we are not concerned by the most likely value of sini > 1 
found by tempo's linear least-squares fit. To confirm that the Shapiro delay is physically 
meaningful, we constructed a map as a function of rric and sini, while optimizing the 
remaining parameters. Our sparse data set is consistent with a large range of allowed pa- 
rameter space, constraining the system to have sini > 0.85 and O.ITMq < tjIc < 0.55 
(95% confidence). The covariance of the poorly-determined rric and sini with the projected 
semimajor axis (asini) and e means that the uncertainties in these parameters given in 
Table 1 are hkely to be underestimated. 

Currently, the strongest constraint on the companion mass is the minimum mass of 
O.195M0 obtained from the mass function by assuming i = 90° and mp = 1.4 M©; the fact 
that Shapiro delay is detected at all suggests that we are likely viewing the system roughly 
edge-on. Future pulsar timing observations near inferior conjunction will significantly im- 
prove the Shapiro delay-derived companion mass. 

2.2. Distance 

We have measured the parallax, vr = 1.22±0.44mas, and hence can estimate the distance 
to PSR J1909— 3744, = 0.82lQ;22kpc. Combined with the proper motion, this distance 
corresponds to a transverse velocity of 140^40 km s"^. -DM-derived distance estimate is 
0.46 kpc (Cordes & Lazio 2002). 

The pulsar's proper motion induces an apparent radial acceleration equal to [j^d/c 
(Shklovskii 1970). This secular acceleration will allow an entirely independent distance 
estimate from precise measurement of the orbital period derivative, Pb (Bell & Bailes 1996). 
For circular binaries such as PSR J1909— 3744, the intrinsic Pb due to general relativity is 
negligible. However, the secular acceleration will eventually produce a measurable Pb which 
can be used to determine the distance. The peak-to-pcak timing signal due to parallax is 
~1.4/is, while the Pb signal increases with the observing baseline, t, as ~ 240 ns (t/1 yr)^. 
Thus, after only a few years of timing Pb will give the more precise distance. We note that 
the acceleration induced by differential Galactic rotation is about two orders of magnitude 
smaller than the secular acceleration. 

The secular acceleration also corrupts the observed period derivative (P). Since the 
pulsar's intrinsic spindown rate (Pint) is non-negative, we can obtain an upper distance limit 
of 1.4 kpc, consistent with the measured parallax. Conversely, the measured proper motion 
and parallax-derived distance allow us to correct P for the secular acceleration and hence 
determine Pint. 
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3. OPTICAL OBSERVATIONS OF THE WHITE DWARF COMPANION 

At the position of the pulsar, the digitized sky survey plates showed a possible coun- 
terpart. To verify this, we obtained images on the night of 2003 June 4 with the Magellan 
Instant Camera (MagIC) at the 6.5-m Clay telescope on Las Campanas. We took a 10-min 
exposure in B and a 5-min exposure in R. The conditions were poor, with thin clouds and 
1" seeing. 

The images were reduced using the Munich image data analysis system (MIDAS), fol- 
lowing the usual steps of bias subtraction (separately for the four amplifiers) and fiat-fielding 
using dome fiats. For the astrometric cahbration, we selected 58 USNO-Bl.O (Monet et al. 
2003) objects that overlapped with our images, were not overexposed, and appeared stellar. 
For these objects, we measured centroids and fitted for zero-point position, plate scale, and 
position angle. Rejecting 23 outliers with residuals larger than O'.'S, the inferred single-star 
measurement error in both bands is O'.'IS in each coordinate, which is consistent with ex- 
pectations for the USNO-Bl.O measurements. Hence our astrometry should be tied to the 
USNO-Bl.O system at the ~0'.'03 level. 

In both images, we found an object near the position of the pulsar which is blue relative 
to other stars in the field; see Fig. 3. The position inferred from our astrometry is aj2ooo = 
19'^09™47!457, 5j2ooo = — 37°44'14'.'17, which is consistent with the timing position listed in 
Table 1 within the uncertainty with which the USNO-Bl.O system is tied to the International 
Celestial Reference Frame (~0'.'2 in each coordinate). 

The identification of this object with the pulsar's companion was confirmed by spectra, 
taken on the night of 2003 June 6, using the Low Dispersion Survey Spectrograph 2 (LDSS2) 
at the Clay telescope. We took two 30-min exposures with the high resolution grism and a 
1" slit, which gives a resolution of ~6 A. The conditions were mediocre, with thin cirrus and 
0'/7 seeing. 

We used MIDAS to reduce the spectra. They were bias-subtracted and fiat-fielded 
using normalized dome flats (in which the noisy blue part was set to unity). After sky 
subtraction, the two spectra were extracted using optimal weighting, and added together. 
To derive the response, we used spectra of Feige 110, taken during the same night, and fluxes 
from calibrated STIS spectra.^ The result is shown in Fig. 4. We stress that since cirrus 
was present, the absolute flux calibration is not reliable, although both from the spectrum 
and from a rough photometric calibration using USNO-Bl.O magnitudes, it follows that the 
companion has F ~ 21. 



^See http://www.stsci.edu/instruments/observatory/cdbs/calspec.html. 
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The spectrum (Fig. 4) shows only the Balmer hnes of Hydrogen, from Ha up to HIO, 
consistent with a DA white dwarf of low surface gravity (and thus low mass). The strength 
of the lines and the slope of the spectrum are similar to what is seen for PSR J1012+5307 
(van Kerkwijk et al. 1996) and PSR J0218+4232 (Bassa et al. 2003), indicating a similar 
temperature of ~8500K.^ 

Assuming a mass of ~O.2OM0, the white dwarf should have a radius of approximately 
0.03 Rq and a luminosity of ~4 x lO'^Lo (Driebe et al. 1998; Rohrmann et al. 2002). The 
corresponding absolute magnitude is My ~ 11.0, which implies V ~ 20.5 at a distance of 
0.8 kpc, consistent with the rough photometry quoted above. 

4. CONCLUSIONS 

Pulse timing precision is often limited by systematic errors that distort the shape of 
pulsar profiles due to the interstellar medium and imperfect polarimetric and instrumental 
response. Intrinsically narrow pulses are less susceptible to these errors, as are pulsars with 
small dispersion measures. PSR J1909— 3744 is likely to become a cornerstone of future 
timing array experiments to detect low-frequency gravitational waves. It is noteworthy that 
PSR J1909— 3744 is far from other high-precision pulsars on the sky and therefore probes a 
unique line of sight. 

In the shorter run, further timing of the pulsar, especially with CPSR2, should yield an 
accurate inclination and a reasonably precise companion mass. Furthermore, in a modest 
spectroscopic observing campaign on an 8-m class telescope, the radial velocity orbit of the 
companion can be measured precisely. These combined results should allow a very reliable 
pulsar mass measurement, and advance our understanding of binary pulsar evolution and 
recycling. 

We thank R. Edwards for invaluable help with pulsar search software, H. Knight for 
assistance with Parkes observations, and M. Hanmy. V. Marino, and M. Roth for help 
with the Magellan observations. BAJ and SRK thank NSF and NASA for supporting their 
research. MHvK acknowledges support by the National Sciences and Engineering Research 
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^This high temperature may be surprising, given the pulsar's long characteristic age of roughly 8 Gyr. It 
indeed sets interesting constraints on the thickness of the outer Hydrogen envelope of the white dwarf; for a 
discussion and references, see Bassa et al. (2003). 
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Table 1. Parameters of the PSR J1909-3744 system. 



Parameter Value'' 



Right ascension, aj2ooo 19''09™47! 44008(2) 

Declination, 5j2ooo -37°44'14'.'226(1) 

Proper motion (mas yr-l) . . . -9.6(2) 

Proper motion in <5, (mas yr~^) .... —35.6(7) 

Annual parallax, w (mas) 1.22(44) 

Pulse period, P (ms) 2.9471080205034(1) 

Period derivative, P (IQ-^O) 1.4026(3) 

Reference epoch (MJD) 52055.8704 

Dispersion measure, DM (pc cm"^) . . . 10.394(1) 

Binary period, Pb (d) 1.5334494503(1) 

Projected semimajor axis, asini (It-s) . 1.897992(4) 

esinuj (xlO-"^) 0.1(90) 

ecosu) (xlO-^) -2.6(17) 

Time of ascending node, T^sc (MJD) . . 52055.87046667(6) 

sini 1-1(4) 

Companion mass ruc (Mq) 0.13(37) 



Derived Parameters 

Mass function, /(m) 

Orbital eccentricity, e (xlO~^) 

Longitude of pcriastron, oj (deg) 

Time of pcriastron. To 

Parallax distance, d^r (kpc) 

Transverse velocity, v± (km s~^) 

Intrinsic period derivative, Pj^t (10"^")'^ 
Surface magnetic field, -Bgurf (xlO^G)*? 

Characteristic age, Tc (Gyr)*". 

Galactic longitude, / (deg) 

Galactic latitude, b (deg) 

Distance from Galactic plane, \z\ (kpc) 

Pulse FWHM, wso (ms) 

Pulse width at 10% peak, wio (/^s) 



''Figures in parenthesis are uncertainties in the last digit quoted. 
Uncertainties are calculated from twice the formal error produced 

by TEMPO. 

''Corrected for secular acceleration based on measured proper 
motion and parallax. 



0.00312196(2) 

2.6(18) 

177.874428 ± 190 
52056.6281373 ± 0.8 
0.82l°;« 

140ifo 

1 q+0.2 
^••^-0.8 
7 7+35.0 
'•'-2.1 

359.73 

-19.60 

0.28l°:i? 

43 

89 
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Fig. 1. — Coherent ly-dcdispersed average pulse profile of PSR J1909— 3744 at 1341 MHz as 
measured by the CPSR2 baseband system. Instrumental smearing is just 2//s. The location 
of a very weak interpulse is indicated by the underscore at left. 




Fig. 2. — Timing residuals for PSR J1909— 3744. (a): Residuals vs. observation epoch, (b) 
Residuals vs. orbital phase. Inferior conjunction occurs at orbital phase 0.16. 
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Fig. 3. — Image of the field of PSR J1909— 3744. The large frame shows our R-band image; 
the position of the pulsar is indicated by the tick marks. The inset shows the B-band image; 
compared to other stars in the field, the companion is blue. 
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Fig. 4. — Spectrum of the companion of PSR J1909— 3744, showing that it is a low-mass 
DA white dwarf. The Balmer hnes from Ha to HIO are indicated. Note that the absolute 
flux calibration is uncertain by about 50%, as cirrus was present during the observations. 
The relative cahbration, however, should be rehable (except longward of ~ 6400 A, where 
the spectrum might be contaminated by second-order light). 



